Harmonic E/B decomposition for CMB polarization maps 
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The full sky cosmic microwave background polarization field can be decomposed into 'electric' 
(E) and 'magnetic' (B) components that are signatures of distinct physical processes. We give a 
general construction that achieves separation of E and B modes on arbitrary sections of the sky at 
the expense of increasing the noise. When E modes are present on all scales the separation of all of 
the B signal is no longer possible: there are inevitably ambiguous modes that cannot be separated. 
We discuss the practicality of performing E/B decomposition on large scales with realistic non- 
symmetric sky-cuts, and show that separation on large scales is possible by retaining only the well 
supported modes. The large scale modes potentially contain a great deal of useful information, and 
E/B separation at the level of the map is essential for clean detection of B without confusion from 
cosmic variance due to the E signal. We give simple matrix manipulations for creating pure E and 
B maps of the large scale signal for general sky cuts. We demonstrate that the method works well 
in a realistic case and give estimates of the performance with data from the Planck satellite. In the 
appendix we discuss the simple analytic case of an azimuthally symmetric cut, and show that exact 
E/B separation is possible on an azimuthally symmetric cut with a finite number of non-intersecting 
circular cuts around foreground sources. 



Observations of fluctuations in the cosmic microwave 
background (CMB) provide information about primor- 
dial inhomogeneities in the universe. One of the most 
interesting questions is whether there was a tensor (grav- 
itational wave) component. CMB polarization measure- 
ments offer a probe of this signal Q, Q, Q via the 'mag- 
netic' component of the polarization map. This offers 
the opportunity to distinguish between different models 
of single field inflation, which generically predict a signif- 
icant amplitude of gravitational waves, and other com- 
peting models which predict nearly zero amplitudes. For 
example the ekpyrotic scenario generically predicts expo- 
nentially small amplitudes of tensor modes, as do infla- 
tion models where density perturbations originate from 
fluctuations in a second scalar field. 

Polarization of the cosmic microwave sky is produced 
by electron scattering, as photons decouple from the pri- 
mordial plasma and during reionization. Gravitational 
waves produce 'magnetic' (B) and 'electric' (E) polar- 
ization components at a comparable level by anisotropic 
redshifting of the energy of photons. Magnetic polariza- 
tion is not produced by linear scalar (density) perturba- 
tions, so detection of a magnetic component would pro- 
vide strong direct evidence for the presence of a primor- 
dial gravitational wave (tensor) component. There is a 
non-linear contribution to the magnetic signal from grav- 
itational lensing of E polarization, though on the scales 
where the tensor mode signal is large the lensing signal 
is sufficiently small that it is negligible for observations 
up to Planck 1 sensitivity. However ultimately the lensing 
signal (and how well it can be subtracted) may provide 
a limit on the amplitude of tensor modes that can be 
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detected U [oj. l| There are also a variety of other pos- 
sible sources for tensor and vector modes generating a B 
mode signature 0, H, 0, EOl , as well as numerous system- 
atics [Tl|. 

Inflationary models generically predict a Gaussian 
spectrum of linear E and B modes, in which case it is 
possible to use the full likelihood function to constrain 
the tensor amplitude without separating into E and B 
modes. However if there are any departures from Gaus- 
sianity, e.g. due to systematics, foregrounds, magnetic 
fields Q , defects 0, |g or unexpected physics, this could 
give misleading answers. So it is useful to have methods 
to separate out pure B modes for robust detections and 
isolation of unexpected features. For small tensor ampli- 
tudes any mixture of E and B modes would be dominated 
by the scalar E signal, and a joint analysis for the ten- 
sor contribution would be only marginally more optimal 
than using only the pure B modes. 

On small scales it is possible to provide excellent sep- 
aration of the E and B mode power spectra using fast 
quadratic estimators 0, These methods provide 

estimators which, when averaged over realizations, give 
zero if there really is no B signal. However in any given 
realization, such as the one we observe, there will in gen- 
eral be a non-zero B signal. This is a potential obstacle to 
clean E/B power spectrum separation on the large scales 
where cosmic variance is most important. 

CMB observations by WMAP have recently indicated 
a substantial optical depth to reionization [14|- If con- 
firmed this would imply that detection of magnetic po- 
larization could be achieved by observation of a relatively 
small number of high signal to noise B modes on large 
scales (corresponding to the horizon size at reionization). 
The large scales are just where cosmic variance is largest, 
and the effect of incompleteness of the map (e.g. due to 
cuts around the galaxy and foreground sources) makes 
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the E and B modes harder to separate. There is there- 
fore clear motivation for methods to cleanly separate the 
E and B modes on the largest scales, avoiding problems 
with cosmic-variance mixing that arises from the use of 
quadratic estimators, and extracting the B modes with- 
out assumptions about their distribution. 

In this paper we discuss various harmonic methods for 
performing E/B separation, equivalent to separation at 
the level of the map. We start in Section [I] with a brief 
summary of the E/B decomposition. In Section UTI we 
review the tensor harmonics and how the E/B mixing 
enters for observations over only part of the sky, and dis- 
cuss the theoretical properties of the coupling matrices. 
In Section IIIII we show how to separate the modes for 
band limited and non-band limited fields with arbitrary 
sky cuts, and show that the method presented in ear- 
lier work |l5( is optimal for general cuts in the all-scales 
limit. For non-band limited fields there are ambiguous 
modes than cannot be separated. We then show how the 
large scale reionization magnetic modes can be extracted 
from a realistic map by retaining only the well supported 
modes. The method is computationally tractable and 
close to exact, and should allow robust detection of the 
large scale magnetic signal with future CMB polarization 
observations. In Section llVl we demonstrate the perfor- 
mance explicitly for a complicated sky cut geometry, and 
provide comparative estimates of the ability of the Planck 
satellite to detect tensor modes using different methods. 
In the Appendix we review some analytic results for az- 
imuthal cuts from Ref. [15| . show that exact separation 
is possible for non-intersecting combinations of azimuthal 
cuts, and give a slightly improved matrix method for ex- 
tracting the pure E and B modes exactly. 

For a method similar to the general method presented 
here, but working explicitly in pixel space, see Ref. [l6j . 
and other related work in Refs. [II El El Q3 . Though 
we only discuss CMB polarization maps explicitly, the 
techniques could be applied to other data, for example 
shear distributions in weak lensing surveys. Indeed B 
modes in the lensing of pre-reionization gas mayprovide 
a good alternative way to detect tensor modes [l9j . 

I. E/B POLARIZATION 

The observable polarization field is described in terms 
of the two Stokes' parameters Q and U with respect to 
a particular choice of axes about each direction on the 
sky. We use spherical polar coordinates, with orthonor- 
mal basis vectors erg and er^. The Stokes' parameters 
define a symmetric and trace-free (STF) rank two linear 
polarization tensor on the sphere 

V ab = \[Q{<r a e <T b e *%<t\) U(v$o$ + a$o* g )]. (1) 

A two dimensional STF tensor can be written as a sum 
of 'gradient' and 'curl' parts 

V ab = V {a V b) PE-e c ia V b) V c P B . (2) 



where V is the covariant derivative on the sphere, an- 
gle brackets denote the STF part on the enclosed in- 
dices, and round brackets denote symmetrization. The 
underlying scalar fields Pe and Pb describe electric and 
magnetic polarization respectively and are clearly non- 
local functions of the Stokes' parameters. One can de- 
fine scalar quantities which are local in the polarization 
by taking two covariant derivatives to form V°V b, Pab = 
(V 2 + 2)V 2 P E and e b c V c V a V ab = (V 2 + 2)V 2 P B which 
depend only on the electric and magnetic polarization 
respectively. For band limited data one could consider 
taking these derivatives of the data to extract the E and 
B components. The problem is that in the neighborhood 
of boundaries it becomes harder to measure the second 
derivatives, so the noise properties are not straightfor- 
ward: taking the derivatives is effectively increasing the 
noise near the boundaries relative to the rest of the map. 

Rather than taking second derivatives it it useful to 
work with integrals over the surface. As an example we 
focus here on the B polarization, since this is of most 
interest. We define the surface integral 

B w = -2 [ dSWe b c V c V a V ab , (3) 

where W is a real window function defined over some 
patch S of the observed portion of the sky. The factor of 
minus two is included to make our definition equivalent 
to that in Ref. [15|. Integrating by parts we have 

B w = V2 f dSW B b *V ab 
Js 

- 2 I dl a (e b a WV c V cb - e b c V c WV ab ) , (4) 

J dS 

where W Ba b = V / 2e c ( a V{ ) ) V C W is an STF tensor window 
function. Thus extraction of the B polarization amounts 
to measuring a well defined surface integral, and two line 
integrals. It is these line integrals that encode the trou- 
blesome aspect of the E/B decomposition. 

II. E/B HARMONICS 

Since the E/B decomposition is inherently non-local, 
it is quite natural to work in harmonic space. The polar- 
ization tensor V a b can be expanded over the whole sky 
in terms of STF tensor harmonics 

Vab = —7= ^2 (Elm Y(f rn)ab + Bi m Yfi m } ab ^j , (5) 
lm 

where Y9\ ab and Yn m \ ab are the gradient and curl tensor 
harmonics of opposite parities defined in Ref. . From 
the orthogonality of the spherical harmonics over the full 
sphere it follows that 

B lm = V2 [ dSY^Tat. (6) 
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In a rotationally-invariant ensemble, the expectation val- 
ues of the harmonic coefficients define the power spec- 
trum: 



( B {lm)' B lm) = Sl'lSm'mPi 



BB 



(7) 



When we have data over a section of the sphere, the ob- 
served data can be encoded in a set of pseudo-multipoles 
Ei m and Bi m obtained by including a window function 
W in the integral of Eq. @J : 



B 



E 

(lm) 1 



)>E(i m y) ,(8) 

where the Hcrmitian coupling matrices are given by 



W, 



(lm) (lm)' 



dS WY^ ab Y^y 



dSWY, 



c* 



(lm)ab (lm 



Y, 



Gab 



(9) 
(10) 



-(lm)(lm)' — 

The matrix W—n m sn m y controls the contamination with 
electric polarization and can be written as a line integral 
around the boundary of the cut defined by W. The ma- 
trices can be evaluated easily numerically in term of the 
harmonic coefficients of the window Wj m using 2 



W. 



±(hmx)(hm2) 



K-irE^V <2 ' 1 + 1)<2 'i;" (2 ' +1) 
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where m = m\ — m%. 

Vectors are now denoted by bold Roman font, e.g. B 
has components B\ rni and matrices are denoted by bold 
italic font, e.g. W± have components W±(i m ^n m y. In- 
cluding the results equivalent to the above for the E po- 
larization, we then have the harmonic relations 



E\ _ / W+ iW-\ /E 
B / ~ {-iW- W+ i IB 



(12) 



Harmonic methods of E/B separation amount to ways of 
solving these equations for linear combinations of the ob- 
served E, B such as to give results that depend only on 
E or B. In principle the matrices are formally infinite, 
though in practice we measure the components of E, B 
up to some finite £ max . However the pseudo- harmonics 
to finite ^ max can still can contain contributions from the 
underlying fields at all £ unless the fields are band lim- 
ited (i.e. there exists some finite ^ max above which all the 
components of E and B are zero), or the coupling ma- 
trices are sufficiently localized. Thus in general the W± 
coupling matrices are rectangular, including the coupling 
to the E and B modes on all scales. 



Our sign conventions follow Ref. Il5l . 



Eigenstructure 

For square coupling matrices the symmetry of the 
spherical harmonics implies that PW P = — Wl for 
the matrix P(i m ){i> m ') = (^l) m *ii'^mm' which satisfies 
P 2 = 1. This implies that eigenvalues of W come in ± 
pairs: 



We = Ae 



W Pe* = -\Pe* 



(13) 



and hence Tr(XV-) = 0. These properties hold for all 
f mal and are a direct result of our window on the sky 
being a real function. For a window function which is 
positive (or zero) everywhere W+ is positive semidefinite, 
and the eigenvalues are bounded between zero and one 
for window functions normalized to lie between zero and 
one. 

In the limit that £ max - > oo the coupling matrix be- 
comes a projection operator onto the section of sky ob- 
served (from now on we assume the window function is 
taken to be zero or one everywhere), so 



W+ iW- 

-iW- W+ 



W+ iW- 
-iW- W+ 



(14) 



This property is ensured by the completeness of the har- 
monics, and shows that 



w W+ + w+w = w 



wi 



Wl = W+. (15) 



For an eigenvector e\ of W+, with W+e\ = \e\, these 
relations imply that 

Wle x = A(l-A)e A 
W+[W-e x ] = (l-A)[W_e A ]. (16) 

It follows that W-&\ is also an eigenvector of W+ with 
eigenvalue 1 — A, and we can define ei_> so that W-e\ = 
\J A(l — A)ei_x- Hence the eigenstructure of W is given 
by 



W-(e x ± ei_x) = ± VA(1 - A)(e A ± ei_ A ). (17) 



Thus two sets of eigenvectors of 
space of W : 



lie in the null 



W+e Q = 
W + ei = ei 



W-e = 
W e l = 0. 



(18) 
(19) 



The set of vectors e! form a basis for the set of sup- 
ported pure-B and pure-E modes and the eo form a ba- 
sis for modes that are zero within the observed regions 
and therefore cannot be measured. The remaining modes 
with A ^ {0, 1} form a basis for a set of ambiguous modes. 

As {1,1'} — > oo the elements of the coupling matrix 
W+(im)(im)' -* W (lm) ( lm y where W is the coupling ma- 
trix for the scalar harmonics (see Ref |2(|). Complete- 
ness of the scalar harmonics implies W is a projection 
matrix and hence has a fraction ~ f s k y unit eigenval- 
ues and ~ (1 — /sky) zero eigenvalues, where / s k y is the 
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fraction of the sky included in the cut. Thus we expect 
a fraction ~ f s ^ y of the eigenvalues of W + to be unity 
and ~ (1 — /sky) to be zero, with the other eigenval- 
ues making up a fraction oc 1 /£ ma . x corresponding to the 
boundary to area ratio. These results only apply in the 
limit that £ m ax — ► oc, however for finite £ m ax ^ 100 we 
expect a large number of eigenvalues of W + close to zero 
or one, corresponding to modes which are either very well 
or very poorly supported over the observed area. 



III. HARMONIC E/B SEPARATION 

To measure the B only we look for a matrix PL such 
that 



(20) 



Assuming such a Pb can be found, we have then have a 
vector Xb of pure-B modes 



B 



W+ 



B (21) 



that have no dependence on E. The projection onto E is 
generated similarly with 



(22) 



where J is the identity matrix. The pure E and B modes 
are then given by 



3) (I) ■ (-) « 



where 



M = Pl 



W+ 
-iW- 



W+ 



(24) 



Band limited case 



Here we assume there exists some 



- for which all 
components of B, and E, are negligible for I > £ max . 
We can perform a singular value decomposition (SVD) 
to write 



w+ 

-iW- 



where 



R 



RUDV^ w RUDV^ 



I 
I 



(25) 



(26) 



(27) 



The matrices U and V are column unitary, and D is di- 
agonal. Since we have observations over only part of the 
sky, many elements of D will be close to zero, indicat- 
ing modes that are not supported over the observed area. 
The tilded variables are constructed by deleting the cor- 
responding columns and rows of U and V', making D a 
smaller square matrix. The approximation can be made 
the same order as the numerical precision by choosing the 
threshold for the diagonal elements of D small enough. 
Premultiplying by tJ* we then have 



DV^ + iHDV^ (28) 



W f E - iHDV^B (29) 



where the Hermitian matrix H is defined by 

iH = CRU. (30) 

The vector DVH contains essentially all the observable 
information about B. Solving we have 



/ RUU i R t 



u^bJ 



I-UU f 



{l-H 2 )DV^B (31) 
(1 - H^DV^E, (32) 



which achieves the E/B separation. It is in the form of a 
projection operator to remove the left-range of the cou- 
pling, followed by a reduction into the basis of partially 
supported modes by multiplication with U' . For band 
limited skies in the absence of noise there is no informa- 
tion loss. However for modes corresponding to non-zero 
eigenvalues of H the signal to noise is decreased rela- 
tive to doing no separation. This is consistent with the 
understanding that the noise on the second derivatives 
of the data needed to perform direct E/B decomposition 
becomes larger as you near the boundary of the region. 

Note that this method is not useful for extracting large 
scale B modes from CMB polarization observations since 
even though B may be effectively band limited at low £, 
there is also expected to be a E signal with a high band 
limit with £ max > 10 3 . Vectors of harmonics are of size 
n = (f max + l) 2 — 4 which makes the method compu- 
tationally infeasible for £ max > 200 without simplifying 
symmetries. On the full sky it is possible to impose a low 
band limit by convolution, however on the cut sky this is 
not possible without mixing information from inside and 
outside the cut. 



Non-band limited case 

We now consider what happens to the above results as 
the band limit is removed. In the limit t max — > oo we can 
use <|15H to show 

H « -D^V^W+W- + W-W+) VD 1 

-» D VD 1 , (33) 
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and from (|25|l that V diagonalizes W+ with VtW+V" = 
D 2 . Taking V to have the related eigenvectors e.\ and 
ei_A in adjacent columns, it also follows from (|1 T|) that 
V^W-V is block diagonal, where the blocks are ei- 
ther zero or off-diagonal 2x2 matrices with eigenvalues 
±^/A(l — A). This implies that H is block diagonal, with 
elements (° J) or zero, and hence that H 2 is diagonal, 

with elements zero or one. Thus each mode in DV^H 
can either be measured exactly (for the zero eigenval- 
ues of H 2 ) or is completely lost by the E/B separation 
(for the unit eigenvalues). The lost modes that cannot 
be separated correspond to the 'ambiguous' modes dis- 
cussed in Ref. ^(| . Since the zero eigenvalues of H 2 are 
determined by the null space of W , in this limit E/B 
separation amounts to projecting out the non-zero eigen- 
values of corresponding to the boundary terms. For 
non-band limited skies in the limit ^ max — > oo the method 
for projecting out the range of W_ advocated in Ref. [l^ 
is therefore optimal for general sky patches, in addition 
to being optimal in the simple azimuthal case analysed 
in detail (see also Appendix 170 . 

In Fig. ^ we show two modes which for a band limit of 
f nra = 300 are pure B, but have non-zero projection into 
the range of W . It is clear that these are dominated by 
a line integral around the boundary, and as such they 
have significantly worse noise than the other modes due 
to the non-zero eigenvalue of H 2 . The line integral is 
sensitive to E power on scales with I > £ max . As the 
band limit increases the eigenvalue of H 2 tends to one, 
and the line integrals can no longer be measured (the 
signal to noise goes to zero). 

For band limited skies one can use knowledge of the 
band limit to use the information from the modes with 
non-zero eigenvalues of W_. However in practice the E 
polarization is effectively non-band limited as far as ex- 
tracting the large scale B signal is concerned. Thus one 
needs to find modes which lie in the null-space W to 
extract pure B modes. Note that it is the non-band lim- 
itedness of E which implies inevitable loss due to mode 
separation, even in the limit of zero noise. 

Extracting low multipoles 

In the case where the E signal is effectively non-band 
limited, but we observe only a finite number n of pseudo- 
harmonics with £ < £ max , a pure B mode e^B can be ex- 
tracted by finding a vector e lying in the left null space of 
Wf° 5 where is in general an n x oo matrix, coupling 
in E on all scales. This requires that 

\JW™\ 2 =et(W+- W? c W+ ot )e = a (34) 

where here the W+ is the Hermitian square n x n ma- 
trix, = (W+, X) and X is n x oo. For a supported 
mode with W+ e = e this criterion is satisfied because 
| e tW+°| > le^W+l and le^Wf ! 2 must be positive or 
zero. Note that W+e = e is sufficient but not necessary 
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FIG. 1: Pure B modes with non-zero projection into the 
range of W-, for band limited E with ^ max = 300. The 
modes are for an azimuthally symmetric patch of sky with 
where # ma x = 20°, and have m = 1 and m — 2. 
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FIG. 2: Typical Ci for a concordance model with a reioniza- 
tion optical depth of r = 0.17. The top line shows the scalar 
j3-mode power spectrum, thick and thin solid lines show ten- 
sor B and E power spectra at an amplitude where tensors con- 
tribute about l/10th of the large scale temperature anisotropy 
[At — 4 x 10 -10 ). The dashed line gives an indication of the 
noise level with Planck, the dotted line show the B spec- 
trum due to gravitational lensing of the scalar E signal. The 
detailed shape of the curves on large scales depends on the 
reionization history. 
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FIG. 3: The functions (thick lines, the coupling to B) 
and Wf (thin lines, the coupling to E) giving the contri- 
butions to the variance for some nearly pure B modes with 
various degrees of support (e = {0.1, 0.01, 0.001}, top to bot- 
tom), constructed with i? max = 30 and the asymmetric sky 
cut discussed in the text. 



for |et | to vanish — in the case of an azimuthal patch 
there is a left null space for finite £ max even though there 
are no vectors satisfying W+e = e (see Appendix^}. 
However in general, without including information up to 
the band limit of the data or special symmetries, this 
cannot be done. The idea here is to perform E/B separa- 
tion for the low multipoles without including data up to 
the band limit (which due to the £^ ax scaling would be 
infeasible with current computers for general patches). 

In general for finite n there will be no fully sup- 
ported modes, however for an eigenvector e of W+ with 
W+e = (l-e)e it follows that l^W™] 2 = 0(e). In other 
words, approximately pure B modes are simply the well 
supported modes. Note that W e w is necessary but 
not sufficient since it does not ensure there is no coupling 
to power from E on scales with I > ^ max in general. 

The signal variance of e^B is given by 



(|e + B| 



BB 



— /~<EE 



(35) 



where 



^^|(etW + ) im | 2 WT=X;|(etW_ 



)hn\ 



(36) 

The scalar contribution to the C EE power spectrum is 
typically between 1-100 times larger than C EB for t < 
100 and levels of B detectable by Planck (sec Fig.^J, and 
remains of the same order of magnitude up to I ~ 2000. 
Taking Cf E ~ const the E contribution to the variance 
is ~ |etW'.f | 2 , thus we require e <C 0.01 for clean sep- 
aration of B. However there is little point removing E 
to levels much lower than the experimental noise, so for 
noise limited observations larger values of e could be used. 

In Fig. we show a few window functions W, and 
W[~ for nearly pure B modes, using ^ max = 30 and the 



realistic cut discussed in the next section. This low value 
of tax = 30 is not sufficient to extract many well sup- 
ported modes (there is only one with e = 0.001), however 
it allows us to compute the rectangular matrices W+ and 
W up to much higher t in order to explicitly show the 
coupling of higher multipoles as a function of the choice of 
e. The expected behavior is demonstrated in the figure, 
with small values of e effectively removing the coupling 
to E on all scales. 



IV. EXAMPLES 

We now demonstrate explicitly the extraction of large 
scale B modes from a realistic non-symmetric sky cut, 
the l kp2' cut 3 used by the WMAP analysis 21]. This 
excludes the galactic region, in addition to a number of 
other foreground sources, in total excluding about 15% 
of the sky, and is shown as the solid area in Fig. This 
may not be a good cut for polarized observations, but 
the cut is realistic in the sense that it does not have any 
artificial symmetries, and therefore is a reasonable test 
case for how B mode extraction can work in practice. 

As shown in Fig. [3 the large scale reionization signal 
with £ < 10 is large and has high signal to noise, at least 
if the optical depth really is > 0.1. We use £ max = 150, 
which is computationally tractable (~ 10GB of memory, 
few days of CPU time) and sufficient to get most of the 
large scale reionization power. The well supported mode 
with e = 0.001 and maximal signal is shown in Fig. 0] 
With more effort (e.g. distributing the computation over 
a cluster or more efficient algorithms) it may be possible 
to push to £ max ~ 300, possibly enough to extract essen- 
tially all of the separable B modes due to primordial grav- 
itational waves. This may be essential if the reionization 
turns out to be less significant, and of course resolving the 
shape of the magnetic power spectrum is a key check that 
the signal really is due to primordial tensor modes. Sepa- 
ration on smaller scales (e.g. as a consistency check, and 
for studying the lensing signal) would have to proceed on 
azimuthally symmetric cuts over clean regions of the sky 
with a number of cuts around foreground sources (see 
Appendix ^) , or rely on quadratic methods that work 
very well for estimating the separated power spectra on 
small scales. 

Identifying the well supported modes is straightfor- 
ward, and requires diagonalization of the Hcrmitian ma- 
trix W+, 



W+ = U+D+U 



(37) 



Since we only need the well supported modes, in practice 
we only need to compute the eigenvectors with eigenval- 



! http://lambda.gsfc.nasa.gov/ 
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FIG. 4: The nearly pure well supported B mode with maximal signal for a r ~ 0.12 reionization model, ^ max = 150, e = 0.001 
and the 'kp2' sky cut discussed in the text. The Mollweide projection is used for the cut, but polarization orientations are not 
projected. 



ues 4 between 1 — e and 1, for some choice of e and £ max . 
For large £ max there are ~ /sky^max such modes, where 
/ s ky is the fraction of the sky included in the window. For 
less large ! max a significant fraction of modes will only 
be partially supported; losing the subset of these modes 
that are not contaminated with E is the price we pay 
for a computationally tractable analysis with manageable 
f mffl . The well supported eigenvectors define a reduced 
column-orthogonal matrix U+ which projects the cut-sky 
harmonics into the space of well supported nearly-pure 
B modes: 

X B = U[b re CJ|B. (38) 

A map of the nearly pure B modes can then be con- 
structed using the harmonics B' = t/+C/}B. Note that 
though computing U+ is time consuming for large £ max , 
once this has been done for a particular cut, construction 
of B maps for different realizations is computationally 
fast. 

In Fig. [5] we show the construction of a B map for a 
simple simulation without noise 5 . The recovered B map 
is not the same as the input map since we are only using 
f max = 150, and some modes are lost due to E/B sepa- 
ration. However it is clear than many of the large scale 



4 We use the LAPACK routine 'ZHEEVR', see 
http : //www . netlib . org/lapack 

5 Further color maps are available at 
http : //cosmologist . inf o/polar/EBsupport .html 



features are present in the recovered map, and that the 
method works usefully well. 



Noise 

In general there will be anisotropic and correlated noise 
on the polarization measurements. The E/B separation 
methods we discussed give E/B clean separation of the 
signal, however any interpretation of the separated modes 
would have to carefully model the actual noise properties 
in the observation under consideration. In general the E 
and B modes will have complicated noise correlations. 
However in the idealized case of isotropic uncorrelated 
noise on the Stokes' parameters the noise properties are 
straightforward |l5j . In particular the well supported E 
and B modes will also have uncorrelated isotropic noise 

(XbX^jv « {UtBE*U + ) N « 0. (39) 

The complications introduced by correlated anisotropic 
noise do not effect the signal E/B separation, which is 
purely determined by the geometry of the sky cut. As 
the noise in a region inside the cut increases, some of the 
modes will become much more noisy. In the limit of infi- 
nite noise in this region, corresponding to an additional 
cut, these modes have zero signal to noise and can be 
removed. This removes combinations of E and B modes 
in such a way that the number of pure B modes that 
can be measured decreases, and new ambiguous modes 
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FIG. 5: Simulation of noise-free well supported B mode extraction (£ max = 150, e = 0.001) for the 'kp2' cut. Top: input B map 
due to tensor modes. Middle: the observed cut sky including E polarization from scalar and tensor modes (At = 4 x 10 -10 ). 
Bottom: the reconstructed well supported nearly pure B modes. Line orientations show the polarization direction, line lengths 
show the magnitude. For each map the lengths are independently normalized. The sky is represented in the Mollweide 
projection (but polarization orientations are not projected), and the beam size is half a degree. 
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(which are combinations of E and B with finite noise) are 
generated. 



Detection by Planck? 

Fig.[|J]shows how the E/B separation method performs 
with regards to the detection probability 6 for magnetic 
polarization due to tensor modes, taking a simple model 
of the Planck satellite as a test case (assuming isotropic 
Gaussian noise, that foregrounds can be subtracted accu- 
rately outside the cut region, and systematics are negligi- 
ble). We assume the primordial tensor modes are Gaus- 
sian with scale invariant power spectrum, with At being 
the variance of the transverse traceless part of the metric 
tensor, and that other cosmological parameters are well 
known. 

The approximate method applied to the asymmetric 
cut does not perform quite as well as an optimal analysis 
with an azimuthally symmetric cut, but the difference is 
not that large — equivalent to reducing a 99% confidence 
detection to about 95%. The asymmetric cut is not ex- 
pected to perform as well because it increases the area ad- 
jacent to a boundary and hence the number of ambiguous 
modes (though the approximate result may be improved 
slightly by taking larger £ roax ). The plot also shows the 
result of applying the non-exact method (retaining just 
the well supported modes) to the azimuthally symmet- 
ric cut, and shows that the results are almost identical 
to performing the exact separation in this case. In the 
exact azimuthal analysis the result is almost completely 
insensitive to £ max , with £ max = 30 results lying on top 
of those shown. This is because of the small number of 
large signal to noise modes on large scales for a signif- 
icant optical depth to reionization. In the asymmetric 
case much larger £ max is needed to obtain well supported 
modes. 

For r > 0.1 it is clear that Planck should have a good 
chance of detecting tensor amplitudes At ^ 10~ 10 , cor- 
responding to an energy scale V 1 / 4 > 10 16 GeV at hori- 
zon crossing during inflation. The sensitivity to tensor 
magnetic modes however depends on the reionization op- 
tical depth, as shown in Fig. [7| If the optical depth 
turns out to be at the lower end, the asymmetric cut 
with € max ~ 150 would not perform so well compared 
to the exact azimuthal case since there would be signif- 
icant power on smaller scales, requiring £ max S> 100 to 
construct all the well supported modes. 



We follow the method used in Ref. Il5l . Due to the high signal 
to noise of a few of the large scale modes this may be somewhat 
suboptimal, and therefore pessimistic. 




Tensor Amplitude At 



FIG. 6: Tensor magnetic polarization detection probabilities 
for Planck at a confidence of 99 per cent (thick lines) and 95 
per cent (thin lines) using nearly pure B modes. Solid lines are 
for the 'kp2' cut using ^ max = 150 and e = 0.001. Dashed lines 
are for exact E/B separation with an azimuthally symmetric 
cut (with the same /sky) and £m&x = 250. Dotted lines are for 
the azimuthally symmetric cut, but using the well supported 
modes method (^ max = 150 and e = 0.001). For all lines we 
assume isotropic Cnoisc = 4.1 x 10~ 4 fj,K 2 and a reionization 
optical depth r ~ 0.12. The vertical dashed line shows the 
tensor amplitude that would contribute about l/10th of the 
large scale temperature anisotropy (approximately the limit 
that can be set without using polarization) . 



V. CONCLUSION 



We have demonstrated explicitly that E/B separation 
is possible on large scales for realistic non-symmetric sky 
cuts. If the reionization optical depth is large almost all 
the detection significance for tensor modes would come 
from the very largest scales, where mode mixing is most 
important. We showed that these modes can be sepa- 
rated in a computationally tractable way by retaining 
only the well supported modes. The ambiguous modes 
that cannot be separated will be dominated by the E sig- 
nal if the noise is low, and whilst useful for analysis of 
the E polarization would not add significantly to a full 
likelihood analysis for magnetic polarization. 

With map-level separation one can also in principle 
study the full distribution of the B signal, for exam- 
ple testing for Gaussianity. Analysing a pure B map 
is essentially the same as analysing a temperature map, 
and much more more straightforward and conceptually 
cleaner than doing a full joint analysis. For noise limited 
observations, when the loss of ambiguous modes loses a 
significant amount of information, an analysis with E/B 
separated variables still provides a valuable cross-check 
on results from a full joint analysis, and may help to 
isolate systematics. 
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Tensor Amplitude At 

FIG. 7: The probabilities for 99 per cent confi- 
dence detections of magnetic polarization with a simple 
model of Planck, for reionization optical depths r = 
{0.2, 0.17, 0.14, 0.11, 0.08, 0.05} (left to right) assuming sharp 
reionization. The dotted line shows the result if there is no 
foreground-distinguishable B mode signal at / < 30. Exact 
B mode separation is used with an azimuthally symmetric 
galactic cut with / s k y = 0.85. 
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APPENDIX A: CIRCULAR BOUNDARIES 

For a circular boundary at constant latitude 9 = 
(i.e. the boundary of an azimuthal patch), modes with 
different m decouple, W_ is block diagonal, and modes 
with different m can be handled separately. The W 
surface boundary integral evaluates to |l5j 



W. 



mdn 



{lm){V m') 



2 m 



[ui(m)u*i(m) + vi(m)vf, (m)] 



where the vectors 



ui{m) = 



vi (m) 



_ /(i-2) 



(; + 2)! ^ S infl- 



Y, 



I in 



'Q-2)! y/8|m|7r(m 2 
{1 + 2)1 sine? 



-Y, 



Irn 



e ■ 

(Al) 



(A2) 



(A3) 



oo limit 



for I > 2 and some arbitrary <f>. In the £ ma 
this is the spectral decomposition (vi(m) and ui(m) are 



7 http : //www. eso . org/science/healpix/ 



normalized and orthogonal), which follows from the re- 
sults: 



sin 2 9 



l=\m\ 
oo 

£ 

l = | m | 



(i + 2)! 



(i + 2y YlmS1A ■ 



8|m|7r(m 2 — 1) 
= 



In i 



d6 V sin 



(A4) 
(A5) 



for \m\ > 2, and 



/— max(2,|m|) 



0-2)! 



2)! 



sin 



d9 



Y, 



lm 



sin f 



(A6) 



for 1 771 1 > 1. In this limit the eigenvalues are —\mj\m\, 
a special case of the general analysis given in section ITT1 
For a single circular boundary and finite ^ max there are 
at most two non-zero eigenvalues per m, with one for 
|m| = 1 and none for m = 0. 

The azimuthal case has the nice property that even 
if E and B are only available for finite t max one can 
project out exactly the cross-contamination on all scales. 
The matrix W truncated at a finite number of rows 
(but retaining an unlimited number of columns) can be 
written as 



W 



U-D-Vl 



(A7) 



where U— is a column unitary matrix with at most two 
columns per m, and hence the range of W_ can be pro- 
jected out using 



I U U [ 



(A8) 



so that 



P E =P_W+E 



P B =P_W+B (A9) 



losing at most two modes per m (per boundary). The ma- 
trix U- can in practice be constructed for each m by nor- 
malization and orthogonalization (e.g. using a SVD) of 
the two vectors ui (m) and vi (m) (which are no longer or- 
thogonal or normalized for finite £ ma . x ). This projection is 
equivalent to the procedure of projecting out the non-zero 
eigenvalues of W- in Ref. ^5|, for \m\ < £ max — 2n, where 
n is the number of boundaries. If W is smaller than it's 
range due to a finite £ max (i-e. for |m| > ^ max — 2n) 
the zero eigenvalues W do not correspond to separated 
modes. 

Two maps consisting of pure E and pure B can be 
constructed simply from 

E' = (I — Z7_L/1)E, B' = (pureE) (A10) 
B' = (J-U_Z7l)B, E'= (pureB). (All) 

If desired a map of the ambiguous modes can also be 
constructed from the remaining modes. The separation 
is computationally trivial in the azimuthal case due to the 
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separability in m, and is possible up to the resolution of 
the experiment (typically £ max > 10 3 ). 

For n co- axial circular boundaries one looses up to 
4n(f max — n) modes due to the E/B separation, and the 
separation remains exact and separable in m. For non- 
co-axial boundaries exact separation is still possible, even 
though modes with different m are mixed due to the ro- 
tation (i.e. with the Wigner-D matrices D l ,). This 
is because the coupling matrix W for the entire sky 
made up of a set of non-intersecting cuts can be written 
as a finite sum of coupling matrices W_ each with finite 
range, and therefore itself has finite range which can be 
projected out exactly. Thus it is still possible to perform 
exact E/B separation at finite € max with azimuthal cuts 
containing a finite number of non-intersecting circular 
cuts around foreground sources. 

To extract a set of E and B modes for likelihood evalu- 
ation one can proceed following Ref. [l^. We briefly re- 
view the method here with one minor enhancement. First 
we diagonalize W+ (taking it to be square) and discard 
the badly supported modes, to write W+ as U+D + U\ 
so 

£>; 1/2 C/{E « I>y 2 [/{E + iD~ 1/2 ulW-B (A12) 
JD; 1/2 C/|B w £>+ /2 t7{B - iD+ 1,2 u\WJE. (A13) 



and construct the matrix U— by deleting the columns of 
U corresponding to non-zero diagonal elements of D . 
It follows from the Hermiticity of XV that U- is made 
up of vectors in the left null-space of the mode-mixing 

~ 1/2 ~V+ 

matrix D + U + XV-. The pure E and B modes are then 
given by 

C/l£>; 1/2 C/|E « 01-D+ /2 [/{e (A15) 
U f _D~ 1/2 ulB « C/lf>+ /2 C/|B. (A16) 

For isotropic uncorrelated noise these modes also have 
isotropic and uncorrelated noise. This method is equiv- 
alent to that presented in Ref. 0] but slightly faster as 

it replaces the SVD of the D + ' U + W- coupling ma- 
trix with a faster diagonalization of a smaller Hermitian 
matrix by using the Hermiticity of XV . Note that as 
mentioned above one should only include modes with 
\m<\ < tax — 2ri, though in practice for galactic cuts 
these modes are badly supported anyway. The coupling 
matrices for an azimuthal cut can be computed efficiently 
using the results given in the appendix of Ref. . 



Then we do the diagonalization 

D+ 1/2 ulW-U + D+ 1/2 = U D U\_ (A14) 
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